Quantum disorder due to singlet formation: The Plaquette lattice 
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I study the order/disorder transition due to singlet formation in a quantum spin system by means 
of exact diagonalization. The systems is build by spin 1/2 on a two-dimensional square lattice 
with two different kinds of antiferromagnetic Heisenberg interactions. The interaction Jp connects 4 
nearest neighbor spins on a plaquette. The interaction J„ connects the plaquettes with each other. 
If Jp = Jn the systems reduces to the simple square lattice case. If one of the interactions becomes 
sufficiently larger then the other the purely quantum effect of singlet formation drives the system 
into a disordered phase with only short range correlations in the plaquettes and a spin gap. I study 
the transition point by evaluating the spin gap and spin-spin correlations. I compare the results 
with previously calculated data from a non-linear a model approach, spin wave theory and series 
expansion calculations. I confirm a critical value of J„ ~ 0.6 for the quantum phase transition point. 
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I. INTRODUCTION 



Low dimensional quantum antiferromagnets have attracted a great deal of attention during the last decade, both 
theoretically and experimentally. Such magnets show a wide variety of magnetic low-temperature behavior, like 
magnetic long range order or spin disorder with or without a spin gap. Especially the question how magnetic 
long range order is influenced by different parameters like frustration, temperature, spin quantum number or lattice 
topology has been addressed. 

A wide variety of different numerical methods like quantum Monte Carlo, renormalization group calculations, series 
expansion or spin wave theory has been developed and used for carrying out research in this field. All these methods 
have their advantages but also some disadvantages. Monte Carlo calculations suffer severe from the sign problem in 
frustrated systems. Series expansions and spin wave calculations may break down at phase transition points. The most 
promising renormalization group theory, the density matrix renormalization group is still limited to one-dimensional 
problems. 

In this respect the exact diagonalization of finite lattices is still a very powerful method because the calculated data 
is a priori exact for the finite lattice under consideration. The exponentially fast growing Hilbert space necessary 
for the calculations limits the enumeration for spin half system to a total of N=36 spins yet. Still this task is 
computationally very demanding and one of the computational challenges of the new century will be the ongoing task 
to push this limit as far as possible by using all the available resources at hand. Exact diagonalization will always be 
one of the forerunners of using all the promised computer power coming up in the future and nevertheless it will be 
not enough to go too far away. 

Despite the mentioned limitations to rather small finite lattices it has been shown that with particular care one 
can draw important conclusions about the ordering behavior of different low-dimensional quantum spin systems, e.g. 
the frustrated ferrimagnet in Id [Q, the kagome lattice in 2d or the Heisenberg antiferromagnet on the 

body-centered cubic lattice in 3d g . 

In this paper I will examine a particular interesting quantum spin system which shows a purely quantum phase 
transition by singlet formation. There is an experimental study of a new compound N a^RbCui{As04)iCl2 under 
way 1^ which suggest that this system might be described by the model I consider here. I also note that a similar 
model has been studied previously with a non-linear a-model approach and spin wave theory Q , plaquette expansion 
[H and Ising series expansion ^. At the end I will compare the results of all this methods with each other. 
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II. THE MODEL 



I consider the Heisenberg antiferromagnet on a two-dimensional square lattice with two different kinds of antifer- 
romagnetic spin exchange parameters Jp > and J„ > 0. 

H ^ Jp"^ SiSj + Jn ^ ksj (1) 
inter intra 

The summation Jointer with Jp takes place between 4 nearest neighbor spins on a particular plaquette. The J2intra 
with J„ takes place between nearest neighbors on different plaquettes (cf. Fig.|l|). The quantum spin Si equals ^. For 
simplicity I put Jp=l for the rest of the paper and take J„ as the changing parameter. 




FIG. 1: The Plaquette lattice with N=16 spins and periodic boundaries. Bold lines denote Jp bonds, thin lines denote Jn 
bonds. 

I note that the Hamiltonian is symmetric if one exchanges J„ and Jp and therefore I consider the range of J„ G [0, 1] 
only. For J„ = 1 this system represents the simple square lattice. 



III. RESULTS 



A. The classical ground state 

The ground state of the classical model where the spins Si can be considered as classical vectors Si is actually 
very simple. Because there is no frustration in the system it will remain in the Neel-state for all values of J„ > 0. 
This classical ground state is illustrated in the left part of Fig.^. Only if J„ = the ground state becomes highly 
degenerated and then the Neel-state is only one of the possible ground states due to the free rotation of the spins in 
one plaquette respectively to any other plaquette. This behavior is depicted in the right part of Fig.^j. 

Of course the question is now, what will happen if quantum fluctuations come into play? I will discuss this in the 
next paragraph by considering the spin gap and spin-spin correlation data for the quantum system on finite lattices 
with N=16,20,32 and 36 spins and periodic boundary conditions. 

Since the Heisenberg Hamiltonian commutes with the square of total spin S^, S = Si, each eigen state of H 
belongs to a certain subspace of the Hilbert space with fixed quantum number S of the total spin (5^ — S{S + 1)). 
To calculate the ground state and the first excited state of H for finite lattices I use the Lanczos technique. 



B. The spin gap A 



First I will analyze the spin gap A 



A — J^mini^min ~^ 1) -^min (•^'mm) ■ (2) 
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FIG. 2: The classical ground state spin configuration of the Plaquette lattice. Left: The Neel-state for J„ > 0, right: J„ — 
with a free rotation of the local axis in each plaquette. 



For bipartite antiferromagnets like the one studied here one can show that the ground state belongs to the 5* = 
subspace and the first excitation is a triplet with 5=1. Considering our model one can immediately note, that for 
Jn = all the plaquettes will be in a singlet state with the total spin of the plaquette being zero. The first excitation 
will then be a state were one plaquette is in its triplet state and a simple analytical calculation then shows: A = 1. 
As stated for J„ = 1 one obtains the simple square lattice case where a Neel-like magnetic long range order without 
a spin gap is well established . For some value of J„ between and 1 I will observe the opening of a spin gap at 
the phase transition point. 

In Fig.^I show the data from exact diagonalization calculations for lattices with up to 36 spins and an extrapolation 
to the thermodynamic limit. 
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FIG. 3: The spin gap A versus Jn- The crosses, squares and circles represent the calculated data for N=16,20, 32 and 36 
sites. The triangles with the errorbars are the extrapolated data with two input sets: N=16-32 and N=16-36. For the dashed 
and dotted lines see text. 

I use a finite size extrapolation A oc N^^ for the entire range of calculated data but will keep in mind that this is 
a reasonable good approximation for J„ > 0.6 only (apparent in the small errorbars) but much less appropriate for 
Jn < 0.5 (visible in larger errorbars). I used two data sets with N=16-32 and N=16-36. They give very similar result 
for Jn > 0.6, but deviate apparently for 0.3 < Jn < 0.5. But even so one can clearly distinguish the two phases, one 
with a finite spin gap and one without and hence I draw two dotted fit lines for J„ < 0.5 and one dashed fit line for 
Jn > 0.6. 

The extrapolated data shows a small finite value for A at J„ = 1 but this is due to the rough approximation 
only. I assume a finite spin gap opens up above any value greater than the vertical dashed line (cf. Fig.^) in the 
thermodynamic limit. The crossing of the two dotted fit lines with the dashed one can then be approximated as a 
value for the phase transition point: Jn'^"^ ~ 0.65 ± 0.05. 
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C. The spin— spin correlation 



In order to further evaluate the phase transition point I study the spin-spin correlation data plotted in Fig.^. Here I 
show the spin-spin correlation between two spins connected by a J„ bond (siSj) intra and between two spins connected 
by a Jp bond {siS j) inter- 
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FIG. 4: The spin-spin correlation (siSj) between two spins connected by a J„ bond (intra) or by a Jp bond (inter) versus J„. 

The intra-plaquette correlation starts from as expected for non-coupled plaquettes and ends at a certain finite 
value for J„ = 1. The inter-plaquette correlation starts from —1/2 (this value can be calculated analytically) and 
ends at the same value as (siSj) intra due to the symmetry at this point. 

Both curves show an inflection point at around J„ w 0.5 which can be assumed to be connected to the phase 
transition. Therefore I carry out a numerical differentiation d{siSj)^^^^^/dJn in order to study this inflection point in 
more detail in the left part of Fig.^. By numerical differentiation this inflection point is transformed into a minimum 
<5mm which is then shown versus the system size in the right part of Figj^. 
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FIG. 5: Left: The differential spin-spin correlation ^'^ versus J„. Right: The minimum position Smin versus the 
inverse system size 1/N. 

In the left part of Fig.^ one can see that with increasing system size the minimum shifts towards greater J„ and 
also gets deeper. One may argue that in the thermodynamic limit one will observe a delta peak right at the quantum 
phase transition point. In the right part of Fig.^ I plot the peak position Smin versus the inverse system size and 
extrapolate a value J!^"^ k, 0.55 ± 0.05. This result is close to the value from the spin gap data. 



IV. CONCLUSION 



In this paper I presented a exact diagonalization study of an antiferromagnetic Heisenberg quantum spin systems 
with s=l/2 with two kinds of interactions, the Plaquette lattice. I calculated the spin gap and the spin-spin correlation 
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for finite lattices and obtained a quantum phase transition at J^'"** « 0.6 ± 0.1. This is a purely quantum phase 
transition due to singlet formation which has no classical counterpart. Our result compares well with 4th order 
plaquette expansion result of Jn'^"* ~ 0.54 [|j or the Ising series expansion result J„'^'"** ~ 0.55 ||^. Our result shows 
clearly that the spin wave result Jn'^"* ~ 0.112 |Q can be ruled out. The non-linear a model approach does 
completely fail to predict any transition to an ordered phase without a spin gap. 

Whether this phase transition is first or second order is still to be answered, the spin-spin correlation data with a 
developing delta peak points more to a first order transition. But clearly a more detailed study is necessary to clarify 
this and corresponding work is under way. 

Anyway, this study has shown that the exact diagonalization of finite lattices can be successfully used for complex 
many-body problems and that the results may help to verify and extend the validity of other methods. 
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